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ABSTRACT 

We investigate the effect of long-range scalar interactions in dark matter (DM) models of 
cosmic structure formation with a particular focus on the formation times of haloes. Utilising 
TV-body simulations with 512 3 DM particles we show that in our models dark matter haloes 
form substantially earlier: tracing objects up to redshift z~6we find that the formation time, 
as characterised by the redshift z 1 / 2 at which the halo has assembled half of its final mass, is 
gradually shifted from z 1 / 2 ~ 1.83 in the fiducial ACDM model to z 1 j 2 ~ 2.54 in the most 
extreme scalar-interaction model. This is accompanied by a shift of the redshift that marks the 
transition between merger and steady accretion epochs from ps 4.32 in the ACDM halos 
to z* ps 6.39 in our strongest interaction model. In other words, the scalar-interacting model 
employed in this work produces more structures at high redshifts, prolonging at the same 
time the steady accretion phases. These effects taken together can help the ACDM model 
to account for a high redshift reionisation as indicated by the WMAP data and can alleviate 
issues related to the survival of the thin-disk dominated galaxies at low redshifts. 
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1 INTRODUCTION 

The A cold dark matter (ACDM) model has proven to be capa- 
ble of explaining a tremendous amount of observational data. As 
we entered the era of precision cosmology — both observationally 
and from a modellers perspective — we are left with a more and 
more detailed picture of what actually took place when the Uni- 
verse cooled down and formed structures starting from z ~ 1100 
to the present day. But there are still issues related to the struc- 
ture formation on galactic scales to be resolved. For in stance, ob- 
ser vations of high red shift quasars dWillott et al.ll20l6T) and galax- 
ies dOuchi et al .1 (2009 ) imply that the epoch of the reionisation of 
the Universe ended at z ~ 6. On the other hand observations 
of t he polarisation of the cosmic microwave background radia- 
tion jKomatsu & et al.l20ld : lBenson et allioOrllKogut et al.l2003h 
suggest that the Universe was already reionised at z;S10 — 11. 
However, galaxy formation within the ACDM model using the sug- 
gested moderate normalisation of the power spectrum of primordial 
density pert urbations crs = 0.8 (as favoured by recent WMAP ob- 
servations; Komatsu & et al. 2010]) may not happen early enough 
to account for such a high redshift of reionisation. Furthermore the 
astoni shing array of recent obser vations of ultra-faint distant ob- 
jects jBouwens et alj|201ol . I2009T) as well as advancements in the 



E-mail: pchela@camk.edu.pl 



understanding of the high-z physics dGonzalez et alj20l^ ; lHaimar] 
120041 : IVolonteri et al.ll2003h have shown that the young Universe 
was a busy and already structure-rich place. In addition, there are 
issues related to the survival of thin-disk dominated galaxies. It 
was shown that in a pure CDM model of galaxy-halo formation, 
the majority of Milky-Way-sized haloes with Mo ~ 10 12 /i _1 M Q 
have accreted at least one object with greater total mass than the 
M ilky Way disk (m > 5 x 10 10 /i -1 M Q ) over the redshifts z < 
3 dDrivedl20T(il; IStewart et al.ll2008h . It ha s been shown recently 
I IGovernato et alj|2009l ; iMoster et alJl201Ch . that the disk may re- 
form after mergers or even survive minor mergers, somewhat eas- 
ing the problem of disk survival. Still, late-time bombardment of 
Milky- Way like haloes, predicted within the ACDM model, could 
still be considered as a concern about the survival of thin-disk- 
dominated galaxies and the existence of galaxies with unusually 
quiet merger history like our Galaxy. Moreover, recent observa- 
tions of the distant Universe suggest that much of the stellar mass 
of bright galaxies was already present in place at z > 1, this 
presents the "downsizing" challenge for models of galaxy forma- 
tion because massive haloes are assembled late in the hierarchical 
clustering process intrinsic to the CDM cosm ology dCowie et al.l 
l2003l : lFontanot et alJl2009l : lBower et alfeOOfj) . Finally In addition 
there were also reports on massive superclusters found at interme- 
diate redshifts (z ~ 1), which seem t o be very rare and extreme 
objects within standard ACDM model tYa mila Yarvura et alj20ld : 
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IColless & et"aill2003l : Iswinbank & et alj|2007t) . Finally there is a 
recent claim of possible observational evidence for a "fifth force" 
in DM sector dLeel2010l) . 

Even though we would not go as far as to call this a 
"galaxy formation crisis", there are noticeable tensions between 
what observations suggest and what (cosmological simulations of) 
the ACDM model seems to predict. The ReBEL (daRk Break- 
ing Equivale nce principLe) model dis cu ssed in recent years in 
the literature dGubser & Peebles! l2004bllat iFarrar & Peeblesll2004 
iFarrar & Rosenl 120071 : IPeebles & Nusseri l2010l) has the potential 
to alleviate these tensions and to help the theoretical ACDM 
model to better match the observations. The ReBEL model is a 
slight modification of the standard cold dark matter model. It in- 
troduces an additional long-range scalar mediated force only in 
the dark sector. This has been studied in numerical simulations 



dNusser et al.l2005l:lHellwing & Juszkiewicd2009l ; lHellwind2oT(il ; 
iKeselman et al.ll201(j| ; lHellwing et alj|201fj) and these studies have 
shown that the model can match observational properties of the 
large-scale structures (like the one-dimensional Ly-a power spec- 
trum and the power spectrum of the SDSsQ galaxies) as well 
as the standard ACDM does, providing at the same time inter- 
esting features with respects to the formation of galactic struc- 
tures. We also acknowledge that recently other authors performed 
stru cture formation tests in similar scalar-interacting DM scenar- 
ios jBaldi et alj|20ld. lBaldill2009l ; iLi & Barrowll201(]| : iLi & Zhaol 
l20ir£ but they focused on models containing scalar-field-like Dark 
Energy interacting with DM. Elementary considerations related to 
the ReBEL model and the reionisation mechanism showed that 
scalar-interacting DM could acco mmodate early reionisation of the 
intergalactic medium JCer]|2006l) . The low-resolution simulations 
presented in the ReBEL literature have indicated that this model 
can foster a picture in which structure formation starts at higher 
redshifts. In this letter we will show that this is indeed the case. We 
will present, for the first time in the literature, results of analysis 
of the high-resolution N-body simulations of hierarchical DM halo 
formation in the ReBEL framework. 



2 SELF-INTERACTING DM: THE REBEL MODEL 

The theo retical ground of the ReBE L model was formed in two 
papers bv lGubser & Peebles ( 2004b a). We follow the phenomeno- 
logical description o f thi s scalar-interacting DM model pr esented in 
iNusser et al] J2005h and lHellwing & Juszkiewiczl J2009h . We con- 
sider a modified effective potential and force law between DM par- 
ticles of the forms 
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where the two free parameters of the model are: /3, the measure 
of the strength of the scalar force compared to usual Newtonian 
gravity, and r s , the screening length of the additional force that is 
constant in a comoving frame. Previous considerations in the liter- 
ature of the subject provide a crude estimate of the values of this 
free parameters to be of orders /3 ~ 1 and r s ~ l/i -1 Mpc. In this 
short report we explore two possible values of the j3 — 0.5, 1 and 
two values of the screening length parameter r s = 0.5, lh^ 1 Mpc. 
We label the corresponding simulations runs as ACDM, b05rs05, 



b05rsl, blrsl , where b stands for /3 and rs marks the r s parameter. 
In the ACDM run we use /3 = (i.e. no scalar forces). Equations 
d 1 12b yield a simple phenomenological description of the ReBEL 
addition to the CDM model. 



3 NUMERICAL MODELLING 

To follow the formation of structures within the ReBE L frame- 
work we use an adapted version of the GAD GET 2 code dSpringell 
120051) . For the detailed descriptions of the modifications made to 
the code we refer the reader to ou r previous paper on this sub- 
ject dHellwing & Juszk iewicz 200%. We have conducted a series 
of high-resolutions DM only A-body simulations containing 512 3 
particles within periodic box of 32ft -1 Mpc comoving width. The 
cosmology used in the simulations was the canonical ACDM with 
f2 m = 0.3, S1a = 0.7, erg = 0.8 and h = 0.7. Thus our mass reso- 
lution is m p ~ 2.033 x 10 7 h^ 1 Mq. The force softening parameter 
was set to be e = 6h~ x kpc. All simulations were started with the 
same initial conditions at Zmi = 50, and varied only in the ReBEL 
parameters. For each run we saved 30 snapshot equally spaced in 
the logarithmic time scale starting from z — 6.092 to z — 0. 

We used the MPI+OpenMP hybrid AMIGA halo finder (AHF) 
to identify haloes and subhaloes i n our simulatioiff . AHF is the suc- 
cessor of the MHF halo finder by iGill et al.l 1 2004 ), a detailed de- 
script ion of AHF is given in the code paper Knollmann & Knebel 
120091) . We note that we have adjusted the code slightly to take the 
modified gravity of the ReBEL models into account. Starting with 
the halo catalogues at z = we construct the merger trees for each 
simulation by cross-correlating the particles constituting the haloes 
in consecutive time steps. For each halo we record all progenitors 
and select the main progenitor as the halo that maximises the merit 



function Ai A^/A^ , where Ni } 2 and N a are the number of parti- 
cles of the haloes and the number of shared particles, respectively. 



4 THE RESULTS 

In this Letter we focus primarily on the mass accretion histories of 
dark matter haloes and their masses, respectively. A more in-depth 
study of the effects of the ReBEL model on their internal properties 
will be presented in a companion paper (Hellwing, Knollman & 
Knebe, in preparation). 



4.1 Mass accretion histories 

Following the scheme presented in IWechsler et al. I d2002h (here- 
after W02), we fit each halo's mass accretion history (MAH) to 
the exponential law of the form 



M(a) = exp[a(l - a )] , M(a) 



M(a) 
M 



(3) 



where M is the halo mass for a given cosmic scale factor a ex- 
pressed in terms of its final mass Mo at z — and a is a free 
parameter to be determined via fitting. 

The fit parameter can be used to express a characteristic red- 
shift of formation: 

«. = --!, (4) 
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AHF is freely available from |http : / / popia ■ ft ■ uam. es/AMIGA| 
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Table 1. The mean and the lcr scatter of the zj and z* distributions. 



Figure 1. The probability distribution of the formation redshift z 1 /2 for all 
haloes containing at least 500 particles (Mq > 10 10 h~ 1 Mq). 



defined as the time when the logarithmic slope of the accretion rate 
dlog Af/dlog a falls below some specified value S. In our studies 
we adopt vale S — 2 advocated in W02. Therefore our 2* is the 
redshift indicating the time when a dark matter halo has entered 
the steady accretion phase of its mass accumulation (i.e. the last 
major merger should have taken place at 2 > The formation 
redshift — usually defined as th e time when an objec t constituted 
half of its present-day mass (e.g. lLacev&Cokll 19931) — can now 
be expressed in terms of 2* 



-log 0.5, -log 0.5 

21/2 = § — z* + 1) = -r — 

b a 



(5) 



We will use 2* and z x /i as obtained by fitting our numerical MAHs 
to eqn. [5] as estimators of the steady accretion transitions and the 
formation times of our halos. We acknowledge that our derivation 
of both these redshifts is obviously not independent. To this extent, 
we also derived the formation redshift Z1/2 by solely using the nu- 
merical MAH of each halo finding the point where the mass drops 
below Mo/2; the cross-correlation between these two Zx/2 values 
clearly follows a straight 1 : 1 relation with a scatter of less than 5%. 
We therefore decided to stick to the value obtained by eqn. (O. 

In general, eqn. l[3} was show n to be a universa l good fit to the 
halos MAH's (e.g. W02), however iTasitsiomi et alj J2004I) pointed 
out that the simple exponential one parameter form of W02 is not a 
good fit to the MAH of cluster haloes as well as haloes that expe- 
rienced recent major mergers. To not be biased by including such 
objects in our study we removed all haloes whose x 2 " va lue devi- 
ated from the mean by more than lcr and all haloes whose Most 
Massive Progenitor (MMP) was not present in at least 24 (out of 
30) snapshots, thus eliminating of order ~ 25 per cent of all haloes 
in the original catalogues. We further limit the analysis to objects 
containing in excess of 500 particles (at redshift 2 = 0) corre- 
sponding to a lower mass cut of ~ 10 10 /i -1 Mq. This leave us 
with ~ 7000 objects in the ACDM sample, and ~ 4700 objects in 
the blrsl ReBELrun. 



4.2 Formation and Characteristic Redshifts 

When plotting the distribution function of formation times Zx/2 in 
Fig. Q] we observe a clear trend for the mean of 21/2 to shift to 
higher redshifts in the ReBEL models. We like to note in passing 
that for the ACDM model our results agree with the ones presented 
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Figure 2. Hierarchical structure formation for all models. The relation be- 
tween the formation redshift (21/2) of objects and their masses M. The 
error-bars represents lcr scatter within the binned distributions. The error- 
bars of the ReBEL model where slightly shifted to the right for clarity. 



in lLin et all d2003l their Fig. 2 and 3), when using the same mass 
bins as in said work. 

We calculated the same kind of distribution for the charac- 
teristic redshift 2* (though not explicitly shown here) and list for 
both distributions the respective means and variances in Tab. Q] 
We note that for the strongest ReBEL model considered in this 
Letter the mean value of zx/2 1S 39 per cent larger then the cor- 
responding A.CDM value. This effect is even more pronounced 
for the mean values of the 2, distributions, the corresponding shift 
is more than 48 per cent there. The corresponding lcr scatters are 
also larger. This is not surprising as these two quantities are linked 
with linear relation of eqn. l[5). We also note that the magnitude of 
changes of (2:1/2) and (2*) in models b05rs05 and b05rsl are of 
comparable order. The reason for this is the following: our sample 
is dominated by objects with masses M < l0 12 h~ 1 M Q whose 
virial radii are r v i r < 200/i~ 1 kpc; therefore for majority of our 
objects the screening length of the ReBEL force is much greater 
than their radii, hence for this class of objects only the (3 parame- 
ter measures overall enhancement of scalar-attraction. To conclude 
this paragraph we denote that for the parameter space probed by 
our simulations the j3 parameter plays a more important role on the 
examined quantities than the screening length. 

Viewing the shifts for formation and characteristic redshift in 
combination we determine that, on average, haloes in the ReBEL 
model form earlier, thus arriving at moderate and low redshifts with 
higher masses. Moreover haloes with exceptionally long quiet ac- 
cretion epochs are much more abundant in the populations of our 
ReBEL models. We will return to the issue of larger masses be- 
low in Sec. 14.31 where we study the halo mass function at various 
redshifts. 

To check whether the hierarchical structure formation scenario 
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Figure 3. The cumulative mass functions for all haloes containing 20 and 
more particles calculated at redshifts z = 6, 3 and 0. Mass functions for 
redshifts z = 6 (z = 3) were scaled down by 10 — 2 (10 — x ) for clarity. 



is preserved in the ReBEL models and how the effect of the earlier 
structure formation depends on the halo mass we show in Fig. [2] 
the average zi/% binned in halo mass. We see that the hierarchical 
character of the structure formation is preserved in all the simula- 
tions runs, i.e. small mass objects form first with higher mass ob- 
jects forming later. In all the ReBEL runs the formation redshifts 
are higher for all considered masses, but the effect is strongest for 
less massive haloes. However, we also notice a clear dependence 
on the parameters of the ReBEL model. For small /3 and small 
r„ the effect of shifting structure formation to the higher redshifts 
is less pronounced. This is of course expected, since the screen- 
ing length acts here as the effective radius of an enhanced accre- 
tion and the /? measure the strength of this enhancement. It is in- 
teresting to note, that for the b05rs05, b05rsl runs the aforemen- 
tioned effect disappears at some fixed mass scale. In other words, 
we see that the averaged ReBEL Zi/ 2 values start to agree with 
the ACDM values starting from mass M ~ 5 x lO^/i^Mo. 
It is only the blrsl model that deviates throughout the whole 
mass range probed by our simulations. Nevertheless we observe 
that in the b05rs05, b05rsl simulations the haloes with masses 
I0 11 h,- 1 M & < M < 5xl0 12 /i _1 M© form at significantly higher 
redshift compared to the standard CDM model. For the blrsl case 
even haloes with masses comparable to a mass of a group of galax- 
ies 10 13 /i -1 Mq < M < 10 14 h~ 1 M e acquire most of their 
mass at relatively earlier times, however we denote that at masses 
M > 1O 13 /i -1 M0 the signal is dominated by small-number ef- 
fects, i.e. we only have of order ~ 30 objects in that mass range: 
to more thoroughly check the differences in ACDM and ReBEL 
models at this mass scale simulations with a larger box size are 
required. 



4.3 The Evolution of the Mass Function 

As speculated before, haloes in the ReBEL models not only form 
earlier but also carry a larger mass. To actually confirm and quan- 
tify this previously made claim (cf. Sec. l4.2t we plot in Fig. [3] the 
cumulative mass functions (CMFs) at the redshifts 2 = 6,3, and 0. 

We notice that for z = 6 and z = 3 all ReBELs CMFs 
have higher amplitudes than the fiducial ACDM CMF in the whole 
probed mass range: At high redshifts the ReBEL model produces 
more non-linear structures compared to the standard CDM model. 
However structure formation in the ACDM paradigm seems to be 



'catching' up with ReBEL at lower redshifts. This can again be seen 
in Fig. [3] the discrepancies of the amplitudes at redshift z — are 
much suppressed in comparison to the high redshift CMF's. In ad- 
dition, we note an interesting feature at the final redshift: while the 
CMFs have a higher amplitude in the ReBEL models for redshifts 
z — 6 and 3 across the whole mass range, at z — only the high 
mass end of the ReBEL mass functions shows an excess relative to 
the fiducial ACDM model. The low mass end of the ReBEL CMFs 
at this redshift show deficiency of haloes compared to the LCDM. 
The latter effect indicates that in the ReBEL model some of the 
low mass haloes were already used in structure formation processes 
as bricks for more massive objects. All this taken together clearly 
points out that the process of rapid non-linear structure formation 
takes place at later times in the ACDM model compared to model 
with the ReBEL forces. 



5 DISCUSSION 

The main conclusion of this Letter is that a slight modification to 
the CDM model in the form of scalar interactions can accommo- 
date more gravitationally bound structures at high redshifts (e.g. 
fig. [3} which in turn promotes an early reionisation of the Universe. 
This can be expressed by examining the r\ factor that describe s the 
efficiency o f reionisation of the Inter Galactic Medium (IGM) JCenl 
I200dl2003h : 



'/ 



f e sc{dfM/dt)e uv 



(6) 



where c* is the star formation factor efficiency, f esc is the ionising 
photon escape fraction, d/M /dt is the halo formation rate, euv is 
the ionising photon production efficiency and C is the gas clumping 
factor. Here the numerator reflects the ionising photon production 
rate while the denominator is the ionising photon destruction rate. 
When r\ exceeds a certain threshold, the IGM becomes ionised. As 
the ReBEL models exhibit a richer structure at z = 6, their for- 
mation rate must have been higher than the one of the ACDM run. 
This would make it easier to reach the threshold of reionisation at 
high reshifts in the ReBEL framework. 

We further have to acknowledge that earlier formation of the 
DM haloes fostered by the ReBEL model, results in shifting a sig- 
nificant part of the violent merger events to higher redshifts. This 
indicates that there are many more haloes with extended steady ac- 
cretion histories than expected in the standard ACDM structure for- 
mation scenarios. This would help thin-disk dominated galaxies to 
preserve their disks at low redshifts and would also imply that the 
exceptionally long period o f peace experienced b y the Milky- Way 
is no more so extraordinary l lHammer et al J 2007b . 

Finally, we have also shown that long-range scalar interactions 
between DM particles increase to some extent the masses of DM 
haloes. That is also an important effect that may help to understand 
the existence of objects like massive superclusters at moderate red- 
shifts. The simulations presented in this Letter however do not en- 
compass a large enough volume to study the formation of these 
massive objects in detail. 

In conclusion, we have shown that the ReBEL model allows to 
alleviate som e of the tension between predictions of the ACDM and 
observations Benson et alj|2006l: iKogut et al. 2003; Driver! 120101 : 
IStewart et all (20081 ; lYamila Yarvura et alj|20ich . however, as we 



only focused on the issue of formation time, it still needs to be in- 
vestigated if and how internal properties of dark matter haloes are 
affected by the ReBEL formalism. One interesting quantity is the 
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halo concent ration parameter, as i t correlated with the halo forma- 
tion time (eg. lWechsler et al .120021) . Preliminary work shows this to 
indeed be the case, with the mean concentration (cw) being 9 (11, 
11, 13) for ACDM (b05rs05, b05rsl, blrsl) for haloes with masses 
M ~ 2 x 10 12 /i _1 Mq, respectively. However the issue of the 
halo concentration parameters require deeper studies and, since a 
more detailed study of the haloes internal properties is in progress, 
we will cover the concentration parameter issue in a forth-coming 
work. 
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